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OPERATING  CHARACTERISTICS  FOR  COMBINER 
WITH  A  DEAD  ZONE  IN  EACH  CHANNEL 


I  INTRODUCTION 

\ 

^ome  data  processing  shortcuts  are  often  required  in  order  to 
keep  the  computational  burden  in  today's  detection  and  tracking 

systems  within  manageable  limits.  One  strategem  to  accomplish 

< 

this  goal  is  to  quantize  the  sigr^Al  levels  at  various  points  in 
the  receiver  processing  chains  Another  is  to  reject  low-level 
quantities,  and  retain  only  the  larger  terms,  in  the  belief  that 
only  the  latter,  will  lead  to  statistically  meaningful  decisions 
on  signal  presence  versus  absence. 

”Here,  we  investigate -oxia^-sJiOh'- technique ,  where  all  levels 
below  a  breakpoint  or  threshold  value  are  rejected,  that  is,  set 
to  zero,  while  those  signal  levels  above  the  breakpoint  are  re¬ 
tained  in  their  full  accuracy^  In  particular,  this  approach  is 
employed  in  each  branch  of  a  combiner,  as  encountered  in  diversi¬ 
ty  or  multiple  ping  transmission.  The  question  to  be  addressed 

.  I 

r~  /~>r  is  the  cost  of  this  data  reduction  procedure,  in  terms  of  the 

additional  signal-to-noj.se  ratio  required  to  maintain  a  desired 
level  of  performance,  as  measured  by  the  false  alarm  and  detec¬ 
tion  probabilities. 
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PROCESSOR  DESCRIPTION 

The  processor  of  interest  is  depicted  in  figure  1.  Received 
inputs  r^,...,rN  are  either  composed  of  noise-only  or  they  all 
contain  signal  plus  noise.  An  example  of  this  situation  is 
afforded  by  a  multiple  ping  transmission,  with  search  on  the 
range  of  a  possible  target.  The  received  signal  in  each  channel 
(if  present)  is  match-filtered  and  square-law  envelope-detected 
at  the  candidate  time  instant  of  suspected  or  hypothesized  peak 
output . 

At  this  point,  instead  of  simply  summing  up  these  multiple 
outputs,  and  in  an  effort  to  reduce  the  amount  of  information 
sent  on  for  further  data  processing,  the  squared  envelope  xn  in 
the  n-th  channel  is  subjected  to  the  nonlinear  operation  depicted 
in  figure  2.  Namely,  all  input  levels  to  the  nonlinearity  below 
breakpoint  (threshold)  value  L  are  replaced  by  zero,  whereas 
those  levels  above  the  breakpoint  are  kept  as  is.  The  breakpoint 
value  L  is  chosen  so  that  a  specified  fraction  F  of  the  input 
data  to  the  nonlinearity  is  passed,  when  noise-alone  is  present 
at  the  inputs;  the  hope  is  that  F  can  be  chosen  very  small,  with¬ 
out  significant  degradation  in  performance.  Finally,  the  output 
of  the  summer  in  figure  1  is  compared  with  output  threshold  T  for 
purposes  of  deciding  on  signal  presence  versus  absence. 
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Figure  1.  Processor  Block  Diagram 


Figure  2.  Nonlinear  Device  Characteristic 
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ANALYSIS  OF  PERFORMANCE 

The  inputs  { r n }  to  figure  1  are  presumed  to  be  statistically 
independent  of  each  other,  whether  signal  is  present  or  not.  The 
squared-envelope  outputs  { xn }  are,  therefore,  also  statistically 
independent  of  each  other,  with  probability  density  function 
Px(u),  which  is  presumed  known  for  both  cases  of  signal  present 
as  well  as  signal  absent.  The  corresponding  cumulative  distri¬ 
bution  function  and  exceedance  distribution  function  are,  respec¬ 
tively, 

u 

Px(u)  =  J  dt  Px(t)  =  Prob( x  <  u), 

—  00 

00 

Qx(u)  -  J  dt  Px(t)  =  Prob( x  >  u).  (1) 

u+ 

The  nonlinear  device  in  figures  1  and  2  is  characterized 
mathematically  by 

(  0  for  x  <  L  ) 

y=<  >  ;  L>0.  (2) 

v.  x  for  x  >  L  J 

Breakpoint  L  is  presumed  nonnegative,  since  the  output  of  the 
squared-envelope  detector  in  figure  1  can  never  be  negative.  The 
fraction  of  data  passed  by  the  nonlinearity  is 

F  -  Prob( x  >  L)  -  Q  (L) .  (3) 

Since  exceedance  distribution  function  Qx  is  known,  this  equation 
can  be  solved  for  the  required  breakpoint  value  L,  once  fraction 
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F  is  specified.  This  calculation  is  done  for  the  noise-only  case, 
since  the  breakpoint  is  desired  to  be  set  for  this  condition. 

Inspection  of  figure  2  immediately  reveals  that  the 
probability  density  function  of  random  variable  yR  is  given  by 

py ( u )  =Px(L)  6 ( u )  +px(u)U(u-L),  (4) 

where  6  and  U  are  the  delta  function  and  the  unit  step  function, 
respectively.  Therefore,  the  characteristic  function  of  random 
variable  yn  is 

00 

fyU)  *  J  du  exp(iSu)  py(u)  « 

—  00 

00 

Px(L)  +  J  du  exp(i^u)  px(u)  .  (5) 

L 

Finally,  using  the  statistical  independence  of  the  system 
inputs,  the  summer  output, 

N 

z  ■  Z  y„  ■  161 

n*=  1 

has  characteristic  function 


r  iN 

r 

r  iN 

fz(U  -  [fyU)j  - 

lpx(L)  + ; 

j  du  exp(Ku)  Px  ( u )  J  . 

(7) 

L 

For  general  given  probability  density  function  px(u),  the  inte¬ 
gral  on  u  in  (7)  can  be  done  efficiently  by  means  of  a  fast  Four¬ 
ier  transform.  Then  the  numerical  evaluation  of  the  exceedance 
distribution  function  of  z,  namely  Q  (u),  can  be  accomplished  by 
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the  techniques  utilized  in  [1,  2,  3].  This  numerical  approach 
would  have  to  be  carried  out  for  both  cases  of  signal  absent  and 
signal  present,  in  order  to  get  the  false  alarm  probability  Pp  as 
well  as  the  detection  probability  PD.  Specifically, 

PF  “  Qz(T;  noise-°nly>' 

Pn  =  Q  ( T ;  signal-plus-noise).  (8) 

L/  Z 

In  essence,  (7)  characteri zes  the  performance  of  the  proces¬ 
sor  in  figures  1  and  2.  The  remaining  effort  is  the  analytical 
and  numerical  manipulation  of  (7)  into  useful  computer  forms  and 
evaluation . 


7/8 
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EXPONENTIAL  EXAMPLE 

STATISTICS  OF  DETECTOR  OUTPUT 

If  the  inputs  { r n }  to  the  processor  of  figure  1  are  Gaussian, 
then  the  squared-envelope  detector  outputs  {xnJ  are  exponentially 
distributed.  We  take  the  probability  density  function  of  xn  to 
be 

f  0  for  u<0  'j 

p  (u)  -  i  l  for  noise-only.  (9) 

x  l  exp(-u)  for  u>0  J 

This  corresponds  to  a  mean  value  of 


J  du  u  Px(u) 


for  noise-only. 


(10) 


This  choice  of  scaling  at  the  detector  output  does  not  constitute 
any  loss  of  generality,  since  absolute  level  obviously  has  no 
effect  upon  the  receiver  operating  characteristics  of  the 
processor  in  figure  1. 

For  Gaussian  signal  also  present  at  the  system  input,  the 
probability  density  function  of  xn  is 

f  0  for  u<0  'j 

p<u)  -  {  >  for  signal-present.  (11) 

l  a  exp(-au)  for  u>0  J 


Here , 


(12) 


where  R  is  the  signal-to-noi se  power  ratio  at  the  matched  filter 
output.  (If  R  ■  0,  then  a  «  1,  and  (11)  reduces  to  (9).)  Thus, 
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any  signal  processing  gains  associated  with  the  filtering  process 
are  incorporated  in  the  value  of  R.  Observe  that  R  is  the 
signal-to-noi se  power  ratio  per  channel  or  per  ping,  not  the 
"total  signal-to-noise  ratio"  at  the  system  output. 

Another  signal  model,  which  also  leads  to  probability  density 
function  (11)  for  the  detector  output,  is  slow  Rayleigh  fading  in 
the  medium  through  which  the  transmitted  pings  traveled.  That 
is,  during  a  single  ping  duration,  the  medium  attenuation  is 
constant,  but  from  ping  to  ping,  the  attenuation  is  statistically 
independent  and  governed  by  a  Rayleigh  probability  density 
function  on  the  received  signal  envelope. 
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CHARACTERISTIC  FUNCTION  OF  OUTPUT  z 


We  will  determine  the  statistics  of  output  z  of  figure  1  for 
the  signal-present  probability  density  function  of  x,  as  given  by 
(11).  The  case  for  noise-only  will  then  follow  immediately  by 
setting  a  »  1. 

The  cumulative  distribution  and  exceedance  distribution 
functions  of  xn  are  given  by  substitution  of  (11)  in  (1),  that  is 


P  (u) 
x 


Qx(u) 


0 

1  -  exp(-au) 
1 

exp( -au ) 


for  u<0 
for  u>0 


for  u<0 
for  u>0 


(13) 


The  characteristic  function  of  random  variable  y  is  obtained  by 
substituting  (11)  and  (13)  in  (5);  thus 


f  (£)  «  1  -  exp(-aL)  +  du  a  exp(i£u-au)  - 

y  j 

L 

-  1  -  B  ♦  B  3  !XP<1>LI  •  ( 14 ) 

where  we  define 

B  *  exp(-aL);  L  >  0.  (15) 


The  characteristic  function  of  output  z  is  given  by  (7)  as 

-  l1  - *  *  B  ]  * 


N 

( 1-B ) N_n  Bn 

n*0 


exp( i ^Ln ) 
(l-H/a)n 


(16) 
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AUXILIARY  FUNCTIONS 


Define  the  set  of  functions 


E(u,n) 


1  for  u  <  0 


CD 

J  dt  t  ^*?*~fc*  “  exp(-u)  en(u)  -  exp(-u)^~Yj 
u  -  _  ..  v  k«0 


for  u  >  0 


(17) 


for  n  >  0.  Here,  we  used  the  partial-exponential  notation  en(u) 
given  in  (4;  6.5.11].  The  expansion  of  the  integral  in  (17)  may 
be  verified  by  repeated  integrations  by  parts. 

Also,  define  the  set  of  normalized  probability  density  func¬ 
tions 


p„<u> 


for  u  <  0  'i 


n-1 


foru>„ 


for  n  >  1, 


PQ(u)  =  6 ( u ) 


for  all  u. 


(18) 


The  corresponding  characteristic  functions  are 


fn(*> 


(1-U) 


n 


for  n  >  0, 


(19) 


while  the  exceedance  distribution  functions  are 


Qn(u)  -  E(u,n-1 ) 

for  all  u,  n  >  1, 

f  1  for  u 

<  0  'I 

Vu)  * 

• 

(20) 

{  0  for  u 

>  0  j 
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EXCEEDANCE  DISTRIBUTION  FUNCTION  OF  OUTPUT  z 

Since  pn(u)  and  fn(£)  are  a  Fourier  transform  pair  for  n  >  0, 
it  follows  that 

- - - —  and  a  p  (au)  (21) 

(l-U/a)n 

are  a  Fourier  transform  pair.  Then  (16)  allows  us  to  determine 
the  probability  density  function  of  output  random  variable  z  as 

N 

pz(u>  *  (l-B)N-n  Bn  a  pn(a(u  -  Ln))  for  all  u,  (22) 

n«0 

where  the  "shift  factor"  u  -  Ln  is  due  to  the  exp(i£Ln)  term. 

This  is  a  useful  expansion,  even  for  large  N,  since  all  the  terms 
are  positive  or  zero;  there  is  no  cancellation,  as  there  would  be 
for  an  alternating  series. 

The  exceedance  distribution  function  of  z  follows  immediate¬ 
ly  from  (22)  as 

N 

°z<u>  "  X]  (n)  <1"B^N_n  B”  Va(u  -  Ln))  for  all  u.  (23) 
n-0 

Again,  this  series  has  no  negative  terms.  Also,  the  Qn  terms  are 
sums  of  positive  quantities,  as  may  be  seen  by  referring  to  (20) 
and  (17).  We  will  be  interested  only  in  u  >  0  in  the  following; 
then  the  n-0  term  in  (23)  is,  by  use  of  (20), 

(1  -  B)N  Q0(au)  -  0  for  u  >  0.  (24) 
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DETECTION  AND  FALSE  ALARM  PROBABILITIES 

We  now  utilize  (8),  (23),  (24),  and  (20)  to  obtain  the  detec¬ 
tion  probability  as 

N 

pD  “  CD  <1"B)N'n  Bn  E(a<T  "  Ln),n-1).  (25) 

n-1 

Here,  B  is  given  by  (15),  and  a  is  given  by  (12). 

The  false  alarm  probability  is  obtained  by  setting  R  -  0, 
that  is,  a  *=  1 : 

N 

Pp  =  ^  (JJ)  (l-F)N-n  Fn  E  ( T  -  Ln,n-1 )  .  (26) 

n=l 

Here,  we  have  utilized  (3)  et  seq.,  (13),  and  the  fact  that  B  in 
(15)  reduces,  for  a  =  1,  to 

exp(-L)  =  Qx(L;  noise-only)  -  F,  (27) 

which  is  the  fraction  of  data  passed  by  the  nonlinearity  in 
figure  1,  for  noise-only.  In  fact,  (27)  allows  us  to  solve 
explicitly  for  the  required  breakpoint  value  L,  for  this  exponen¬ 
tial  example,  as 

L  -  -ln( F ) .  (28) 

To  summarize,  (25)  and  (26)  give  the  detection  and  false 
alarm  probabilities  in  terms  of  fundamental  quantities 
N,  number  of  channels, 

F,  fraction  of  data  passed, 

R,  signal-to-noi se  power  ratio  per  channel, 

T,  output  threshold.  (29) 
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The  remaining  variables  in  (25)  and  (26)  are  given  by  (28),  (12), 
and  (15)  as 

L  *  -ln(F),  a  -  ,  B  -  exp(-aL),  (30) 

in  order. 


SPECIAL  CASES 


For  N  «  1 ,  one  channel,  (25)  and  (26)  reduce  to 


PF  =  F  Qx ( T-L ) 


F  for  0  <  T  <  L 

exp(-T)  for  L  <  T 


PD  -  B  Q1(a(T-L) ) 


Fa  for  0  <  T  <  L 

exp(-aT)  for  L  <  T 


(31) 


That  is,  PD  «  Pfa  for  N  *  1,  independent  of  the  value  of  fraction 
F.  This  is  obvious  from  (13)  in  this  case. 

Instead,  if  fraction  F  -  1,  that  is,  no  nonlinearity,  then 
(28)  and  (15)  yield  L  -  0 ,  B  *=  1 ,  and  we  find 

PF  -  E ( T , N— 1 )  \ 


\  for  F  -  1.  (32) 


PD  -  E(aT,N-l)  ) 


These  results  agree  with  (5;  (7)  and  (8)).  A  program  for  the 
evaluation  of  general  results  (25)  and  (26),  as  well  as  the 
special  case  (32),  is  presented  in  the  appendix. 
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GRAPHICAL  RESULTS 

* 

In  figure  3  ,  the  receiver  operating  characteristic  (ROC)  is 
given  for  N  *  1,  F  «  1.  That  is,  there  is  one  channel  and  the 
nonlinearity  is  not  active.  There  is  no  need  to  consider  values 
of  F  less  than  1,  according  to  the  comment  under  (31);  however, 
see  the  subsection  below  on  achievable  false  alarm  values.  The 
curves  in  figure  3  are  parameterized  according  to 

R(dB)  «=  10  log  R.  (33) 

The  remaining  fundamental  quantity,  threshold  T  in  (29),  has  been 
eliminated,  and  PD  is  plotted  versus  Pf  on  normal  probability 
paper . 

In  figures  4,  5,  6,  7,  the  number  of  channels  is  kept  at 
N  =  2,  while  fraction 


F  -  1,  .1,  .01,  .001,  (34) 

respectively.  Additional  cases  for 

N  -  4,  6,  8,  16,  32,  64,  (35) 

in  figures  8  through  31,  complete  the  coverage  in  a  similar 
fashion. 


♦Figures  3  through  31  are  grouped  at  the  end  of  this  section. 
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ACHIEVABLE  FALSE  ALARM  VALUES 

Not  all  values  of  false  alarm  probability  can  be  reached  by 

the  processing  system  of  figure  1.  Since  the  nonlinear  device 

output  yn  can  only  take  on  the  values  yn  -  0  and  yn  >  L,  the  sum 

z  can  only  assume  the  values  z  -  0  and  z  >  L.  Also,  since  the 

N 

probability  of  z  «  0  is  (1  -  F)  for  noise-only,  where  F  is  the 
fraction  of  data  passed  by  the  nonlinearity  in  each  channel,  then 

M 

the  probability  of  getting  z  >  L  is  1  -  (1  -  F)  .  Thus,  the 
range  of  reachable  false  alarm  probabilities  is 

PF  <  1  -  (1  -  F)N.  (36) 

This  bound  holds  regardless  of  the  form  of  the  probability  densi¬ 
ty  for  random  variables  (xn)  in  figure  1. 

For  the  special  case  of  N  ■  1,  this  rule  yields  PF  <  F.  Thus, 
the  plot  in  figure  3  for  N  «  1,  F  =  1  must  be  modified  for  F  <  1, 
to  the  extent  that  only  the  values  for  Pp  <  F  are  achievable. 

For  N  >  1,  the  rule  in  (36)  first  becomes  obvious  in  figure  6 
for  N  *=  2 ,  F  «=  .01.  Namely,  (36)  yields 

PF  <  1  -  (1  -  .01)2  -  .0199.  (37) 

Thus,  the  curves  in  figure  6  are  terminated  to  the  right  of  this 
value  of  the  false  alarm  probability.  This  termination  feature 
occurs  in  numerous  other  figures,  always  governed  by  (36). 
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ERRATIC  BEHAVIOR  OF  RECEIVER  OPERATING  CHARACTERISTICS 

Some  of  the  curves  develop  significant  kinks  for  larger 
values  of  the  false  alarm  probability;  see  figure  31  for  the  most 
pronounced  example  in  this  set  of  results.  This  behavior  is  not 
due  to  computer  round-off  error;  rather,  it  is  due  to  the  shifted 
components  of  probability  density  function  (22)  "kicking  in"  when 
the  output  threshold  reaches  various  multiples  of  breakpoint  L. 
Equivalently,  the  shifted  E-function  components  of  the  detection 
probability  and  false  alarm  probability  in  (25)  and  (26)  are 
activated  at  different  threshold  levels,  reflecting  the  inherent 
abrupt  change  of  behavior  of  these  functions  at  zero  argument. 

For  example. 


E(u,0) 


E(u,l) 


1 

exp( -u ) 

V 

( 

1 

exp( -u ) ( 1+u ) 


for  u  <  0 
for  u  >  0 

4 

for  u  <  0 
for  u  >  0 


(38) 


Thus,  E(u,0)  has  a  discontinuous  slope  at  u  «  0,  while  E(u,l)  has 
a  discontinuous  second  derivative  at  u  ■  0. 
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OBSERVATIONS 

The  required  signal-to-noi se  ratios  for  various  values  of  N 
and  F  are  presented  in  tables  1  and  2  for  two  different  levels  of 
performance,  as  read  directly  from  figures  3  through  31.  The 
overriding  impression  is  that  the  degradation  in  performance  is 
not  severe,  even  for  small  values  of  F,  the  fraction  of  data 
passed  by  the  nonlinearity.  For  example,  from  table  1,  the 
decibel  difference  at  F  ■  .001  versus  F  ■>  1  is,  for  N  -  1 ,  2,  4, 
6,  8,  16,  32,  64,  respectively,  just 

0,  0.3,  0.8,  1.0,  1.2,  1.7,  2.2,  2.7  dB. 

For  table  2,  these  differences  are  substantially  the  same: 

0,  0.2,  0.5,  0.8,  1.0,  1.5,  2.0,  2.5  dB. 

Thus,  the  losses  increase  from  0  dB  at  N  -  1  channel,  to  less 
than  3  dB  for  N  =  64  channels. 

The  situation  is  slightly  worse  for  the  lower-quality  case  of 
Pp  =  IE-3,  PD  =  .5.  Namely,  as  F  is  changed  from  1  to  .001,  the 
required  increment  in  signal-to-noi se  ratio  is  0  dB  for  N  ■  1, 
whereas  it  is  3.4  dB  for  N  -  64. 
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N 

1 

Required  R(dB)  for  F  - 
.1  .01 

.001 

1 

12.8 

12.8 

12.8 

12.8 

2 

9.5 

9.6 

9.7 

9.8 

4 

6.8 

6.9 

7.2 

7.6 

6 

5.4 

5.5 

5.9 

6.4 

8 

4.5 

4.6 

5.1 

5.7 

16 

2.3 

2.6 

3.2 

4.0 

32 

0.4 

0.7 

1.6 

2.6 

64 

-1.4 

-1.0 

0.1 

1.3 

Table  1.  Required  Signal-to-Noi se  Ratio  for  Pp  -  IE-6,  PD  *  . 5 


N 

1 

Required 

.1 

R(dB)  for  F  - 
.01 

.001 

1 

22.4 

22.4 

22.4 

22.4 

2 

16.0 

16.0 

16.1 

16.2 

4 

11.6 

11.6 

11.8 

12.1 

6 

9.4 

9.6 

9.9 

10.2 

8 

8.1 

8.2 

8.6 

9.1 

16 

5.3 

5.6 

6.1 

6.8 

32 

2.9 

3 . 3 

4.0 

4.9 

64 

0.8 

1.2 

2.2 

3.3 

Table  2.  Required  Signal-to-Noise  Ratio  for  P^  -  lE-8,  PD  *  .9 
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Probability  of  False  Rlarm 
Figure  3.  ROC  for  N=l,  F=1 
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.0 


Probability  of  False  Rlarm 
Figure  4.  ROC  for  N=2 ,  F=1 


Probability  of  False  Alarm 
Figure  5.  ROC  for  N=2,  F=. 1 
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Probability  of  False  Alarm 
Figure  7.  ROC  for  N=2,  F=.001 
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Probability  of  False  Rlarm 
Figure  8.  ROC  for  N=4 ,  F=l. 
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i  VI  - 

E-10E-9  E-8  E -?  E-G  E-5  E-4  E-3  .01.02  .05  . 


Probability  of  False  Rlarm 
Figure  9.  ROC  for  N=4 ,  F=.l 
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E-10E-9  E-0  E-7  E-S  E-5  E-4  E-3  .01.02  .05 

Probability  of  False  Rlarm 
Figure  io.  ROC  for  N=4 ,  F=.01 
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Probability  of  False  fllarm 
Figure  12.  ROC  for  N=6,  F-\ . 

31 
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Probability  of  False  Alarm 
Figure  14.  ROC  for  N=6,  F“.01 


Probability  of  Detection 


TR  8595 


E-10E-9  E-B  E-7  E-6  E-5  E-4  E-3  .01.02  .05  .1 

Probability  of  False  fllarm 
Figure  15.  ROC  for  N=6,  F=.001 
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Probability  of  False  Rlarm 
Figure  16.  ROC  for  N=8,  F=l. 
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Probability  of  False  Alarm 
Figure  17.  ROC  for  N=8,  F=. 1 
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Probability  of  False  Alarm 
Figure  18.  ROC  for  N=8,  F=.01 
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Probability  of  False  Rlarm 
Figure  22.  ROC  for  N=16,  F«.01 
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Probability  of  False  Alarm 
Figure  23.  ROC  for  N=16,  F=.001 
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*  E-10E-9  E-8  E-7  E-6  E-5  E-4  E-3  .01.02  .05  . 


Probability  of  False  Alarm 
Figure  24.  ROC  for  N=32,  F=l. 
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•  VI  - 

E-10E-9  E-8  E-7  E-6  E-5  E-4  E-3  .01.02  .05  . 

Probability  of  False  Rlarm 
Figure  26.  SQC  for  N=32,  F=.01 
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’  E-10E-9  E-8  E-7  E-6  E-5  E-4  E-3  .01.02  .05  . 

Probability  of  False  Rlarm 
Figure  27.  ROC  for  N=32,  F=.00 1 


Probability  of  False  Rlarm 
Figure  28.  ROC  for  N=64 ,  F=l. 


Probability  of  Detection 


TR  8595 


E-10E-9  E-8  E-7  E-B  E-5  E-4  E-3  .01.02  .05  . 


Probability  of  False  Alarm 
Figure  29.  ROC  for  N=64,  F=.l 
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Probability  of  False  Alarm 
Figure  30.  ROC  for  N=64,  F=.01 
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Probability  of  False  Rlarm 
Figure  3i.  ROC  for  N®64,  F=.001 
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SUMMARY 

The  cost  of  suppressing  the  low-level  outputs  of  the  detected 
squared-envelopes  is  generally  minimal,  unless  the  number  of 
channels  becomes  very  large.  This  conclusion  has  been  drawn  only 
for  the  example  where  these  squared-envelopes  have  an  exponential 
probability  density  function  for  both  the  noise-only  as  well  as 
the  signal-plus-noise  cases.  It  should  also  be  checked  out  for 
other  candidate  forms  of  probability  density  functions  besides 
exponential . 

One  line  of  reasoning  that  makes  this  conclusion  more 
acceptable  is  that  it  is  only  the  larger  outputs  from  the 
detectors  that  are  going  to  lead  to  positive  statements  about 
signal  presence.  Thus,  suppression  of  the  smaller  outputs  should 
be  inconsequential,  at  least  for  few  channels.  However,  for  a 
large  number  of  channels,  the  sum  of  many  nonzero  low-level 
quantities  may  add  up  to  a  significant  value  and  affect  an 
occasional  detection  decision. 
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APPENDIX  A.  PROGRAM  FOR  RECEIVER  OPERATING  CHARACTERISTICS 


10 

!  GENERATE  PD  VS  PF;  COMBINER  WITH  DEAD  ZONE  IN  EACH  CHANNEL, 

TR8595 

20 

N=64  ! 

NUMBER  OF  CHANNELS;  N>=1 

30 

F=.001  ! 

FRACTION  OF  DATA  PASSED;  0.<F<=1. 

40 

DIM  T <100)  ! 

THRESHOLD  VALUES 

50 

COM  Pf <100),Pdl<100),Pd2<100),Pd3<100),Pd4<100),Pd5<100> 

60 

COM  Pd6< 108) , Pd7< 100) ,Pd8< 100) , Pd9< 100) , Pdl0< 100) , Pdl 1 < 100) 

70 

COM  Pdl2<100),Pdl3<100>, 

Pdl4<100),Pdl5<100>,Pdl6(100),Pdl7(100) 

80 

COM  Pdl8<100),Pdl9<100), 

Pd20< 100) 

90 

DOUBLE  N, I, J  ! 

INTEGERS 

100 

T=.  01 

110 

T  =  T  + .  0 1 

120 

Pf *FNPf <  T , F ,  N ) 

130 

IF  Pf >. 1  THEN  110 

140 

T 1 =  MAX  <  T- . 0 1 , .01) 

150 

T=T+ . 0 1 

160 

Pf=FNPf <T,F,N> 

170 

IF  Pf>lE-10  THEN  150 

180 

T2  =  T 

190 

De1t=<T2-Tl )/100. 

200 

FOR  1=0  TO  100 

210 

T  =  T  1  +De  1 1  •*  I 

220 

T  < I >  =  T  ! 

THRESHOLD  VALUES 

230 

Pf  < I ) =FNPf <  T , F , N  )  ! 

FALSE  ALARM  PROBABILITIES 

240 

NEXT  I 

250 

Rldb=-5  ! 

STARTING  SIGH AL-TO- NOISE  RATIO  <dB) 

260 

Delr=.5  ! 

INCREMENT  IN  SNR  <dB) 

270 

FOR  J=1  TO  20 

280 

Rdb=Rldb+<J-l )*De1 r 

SIGN AL-TO- NOISE  RATIO  FER  CHANNEL 

<  dB  ) 

290 

R=10. M . l*Rdb)  ! 

POWER  RATIO 

300 

FOR  1=0  TO  100 

310 

T  =  T(  I  > 

320 

Pd=FNPd<R, T,F,N) 

DETECTION  PROBABILITIES 

330 

IF  J=1  THEN  Pdl < I )=Pd 

340 

IF  J =2  THEN  Pd2< I )=Pd 

350 

IF  J =3  THEN  Pd3< I )=Pd 

360 

IF  J  =  4  THEN  Pd4 < I ) =Pd 

370 

IF  J =5  THEN  Pd5< I )=Pd 

380 

IF  J=6  THEN  Pd6< I )=Pd 

390 

IF  J  =  7  THEN  Pd7 < I ) =Pd 

400 

IF  J=8  THEN  Pd8< I )=Pd 

410 

IF  J=9  THEN  Pd9< I )=Pd 

420 

IF  J=10  THEN  Pdl0< I )=Pd 

430 

IF  J  =  ll  THEN  Pd  1 1 < I ) =Pd 

440 

IF  J=12  THEN  Pdl2< I )=Pd 

450 

IF  J= 1 3  THEN  Pdl3< I )=Pd 

460 

IF  J= 1 4  THEN  Pdl4< I )=Pd 

470 

IF  J= 1 5  THEN  Pdl5< I >=Pd 

480 

IF  J=16  THEN  Pdl6< I )=Pd 

490 

IF  J=1 7  THEN  Pdl7< I )*Pd 

500 

IF  J= 1 8  THEN  Pdl8< I )=Pd 

510 

IF  J= 1 9  THEN  Pdl9< I )=Pd 

520 

IF  J=20  THEN  Pd20< I )*Pd 

530 

NEXT  I 

540 

NEXT  J 
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550  FOR  1=0  TO  100 

560  Pf  <  I  >=FNlrwphi  CPf  <  I  )  > 

570  Pdl<l)=FNInvphi <Pdl(I)) 

580  Pd2<I>=FNIrwphi <Pd2<I>> 

590  Pd3< I )=FNInvphi <Pd3<I>> 

600  Pd4< I )=FNInvphi <Pd4<I>> 

610  Pd5< I )=FNInwphi <Pd5< I > > 

620  Pd6< I )=FNInvphi <Pd6<I>> 

630  Pd7< I >=FNInwphi <Pd7< I ) > 

640  Pd8< I )=FNInwphi (Pd8( I ) ) 

650  Pd9<  I  >  =  FNIr»wphi  <Pd9< I >  > 

660  Pdl0CI)=FNlnvphi <Pdl0<I>) 

670  Pdl 1 < I >=FNInophi <Pdl 1 <  I )  > 

680  Pdl2< 1 )=FNInwphi <Pdl2<I>> 

690  Pdl3< I )=FNInwphi <Pdl3<I>) 

700  Pdl4< I )=FNInwphi <Pdl4< I > > 

710  Pdl5< I >»FNInvphi <Pdl5(I>> 

720  Pdl6< I >=FNInwphi <Pdl6<I>> 

730  Pdl7< I >=FNInophi <Pdl7<I>) 

740  Pdl8< I >=FNInvphi <Pdl8<l>> 

750  Pd  1 9 < I >=FN I nvph i < Pdl9 < I >  > 

760  Pd20< I )=FNInvphi <Pd20( I ) ) 

770  NEXT  I 

780  CALL  Plot 

790  END 

800  ! 

810  DEF  FNInvphi <X>  !  RMS  55,  26.2.23 

820  IF  X= . 5  THEN  RETURN  0. 

830  P=MIN<X, l.-X> 

840  T=-LOG  <  P ) 

850  T =SQR  <  T  +  T  > 

860  P=l.+T*<1.432788+T*<. 1 89269+T» . 00 1 308 > > 

870  P=T-<2.515517+T*<.802853+T*.010328>>/P 

880  IF  XC.5  THEN  P=-P 

890  RETURN  P 

900  FNEND 

910  ! 

920  DEF  FNECU, DOUBLE  N>  !  N>=0 

930  DOUBLE  K  !  INTEGER 

940  IF  U<=0.  THEN  RETURN  1. 

950  S=T=EXP<-U) 

960  IF  N=0  THEN  RETURN  S 

970  FOR  K=1  TO  N 

980  T=T*U/K 

990  S=S+T 

1000  NEXT  K 

1010  RETURN  S 

1020  FNEND 

1030  ! 
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1040  DEF  FNPf<T, F,  DOUBLE  N )  !  FRLSE  RLRRM  PROB.  T >=0 , 0< F<  =  1 , N >= 1 
1050  DOUBLE  Ns,Nl  I  INTEGERS 

1060  IF  F< 1 .  THEN  1090 

1070  Pf =FNE<T,N-1 > 

1080  RETURN  Pf 

1090  L  =  -LOG  <  F  > 

1100  N 1 =N+ 1 

1110  F 1 = 1 . -F 

1120  fl=F/Fl 

1130  T  n = F 1  N 

1140  C=T 

1150  Pf  =  0. 

1160  FOR  Ns* 1  TO  N 

1170  C=C-L 

1180  Tn=Tn*fi*<Nl-Ns>/Ns 

1190  Pf=Pf+Tn*FNE<C,Ns-l) 

1200  NEXT  Ns  • 

1210  RETURN  Pf 

1220  FNEND 

1230  ! 

1240  DEF  FNPd<R,T,F, DOUBLE  N>  !  DETECTION  PROB.  R >=0 , T >=0 , 0< F< = 1 , N >= 1 

1250  DOUBLE  Ns,Nl  !  INTEGERS 

1260  fls= 1 . / < 1 . +R )  !  a 

1270  IF  F< 1 .  THEN  1300 

1280  Pd  =  FNE  <Rs*T, N-l > 

1290  RETURN  Pd 

1300  L*-LOG<F> 

1310  N 1  *  N  + 1 

1320  B=EXP  < -fls*L) 

1330  B 1 = 1 . -B 

1340  fl=B/Bl 

1350  Tn=Bl^N 

1360  C=T 

1370  Pd=0. 

1380  FOR  Ns= 1  TO  N 

1390  C=C-L 

1400  Tn=Tn*fl*<Nl-Ns)/Ns 

1410  Pd=Pd+Tn*FNE<fis»C,Ns-l > 

1420  NEXT  Ns 

1430  RETURN  Pd 

1440  FNEND 

1450  ! 
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1 460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
I860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 


SUB  Plot  !  PLOT  PD  VS  PF  ON  NORNRL  PROBABILITY  PRPER 
COfl  Pf  <  * ),  Pd  1  <«>  ,  Pd2  <  *  > , Pd3  < «  > ,  Pd4  (*>,  Pd5  (*> 

CON  Pd6<*>,Pd7<*>,Pd8<*),Pd9<*>,Pdl0<*>,Pdl 1<*> 

CON  Pdl2(*>,Pdl3<*>,Pdl4<*>,Pdl5<*>,Pdl6<*>,Pdl7<*> 

CON  Pdl8<*),Pdl9<*>,Pd20<*> 

Din  A$£30],B$C30],C$C31I 

Din  Xlabe1*<l j 30> , Y1 abel *< 1 j 30> 

Din  Xcoord<ls30),Ycoord<l!30> 

Din  Xgrid<l:30>, Ygrid<li30> 

DOUBLE  N, Lx, Ly , Nx, Ny, I  !  INTEGERS 
! 

A$="Probab i 1 i ty  of  False  Alarm" 

Bt="Probabi 1 i ty  of  Detection" 

C$="F i gure  ROC  for  N=64,  F=.001“ 

! 

Lx  =  1 2 

REDin  Xlabel$<l:Lx),Xcoord<l:Lx> 

DATA  E-10,E-9,E-8,E-7,E-6,E-5,E-4,E-3, . 0 1 , . 02 , . 05 , . 1 
READ  XI abel $<*) 

DATA  IE-10,  IE-9,  IE-8,  IE-7,  IE-6,  IE-5,  IE-4, . 001 , . 01 , . 02, . 05, . 
READ  XcoordC*) 


Ly=  1 8 

REDIN  Ylabelt(l!Ly),Ycoord<l:Ly> 

DATA  . 01, . 02, . 05, .1,. 2, .3, .4, .5, .6, .7, .8, .9 
DATA  .95, .98, .99, .995, .998, .999 
READ  Y 1 abe 1 ♦ <  * ) 

DATA  . 01, . 02, . 05, .1,. 2, .3, .4, .5, .6, .7, .8, .9 
DATA  .95, .98, .99, .995, .998, .999 
READ  Ycoord<* ) 

! 

N  x  =  1 4 

REDIN  Xgrid(ltNx) 

DATA  IE- 10, IE-9,  IE-8,  IE-7,  IE-6,  IE-5,  IE-4 
DATA  . 001 , . 002, . 005, . 01 , . 02, . 05, . 1 
READ  Xgr i d<»> 

! 

Hy=  1 8 

REDin  Ygrid(liNy) 

DATA  .01, .02, .05, . 1, .2, .3, .4, .5, .6, .7, .8, .9 
DATA  .95, .98, .99, .995, .998, .999 
READ  Ygri d<*> 

! 

FOR  1=1  TO  Lx 

Xcoord< I )=FNInvphi (Xcoordt I > ) 

NEXT  I 

FOR  1=1  TO  Ly 

Ycoord< I )=FNInvphi <Ycoord<I>) 

NEXT  I 

FOR  1=1  TO  Nx 

Xgr i d< I >=FNInoph i <Xgr i d< I ) ) 

NEXT  I 

FOR  1=1  TO  Ny 

Ygr i d< I >  =FNInophi < Ygr i d< I >  > 

NEXT  I 


56 


TR  8595 


2010  Xl=Xgr i d<  1 ) 

2020  X2  =  Xgr i d  <  Nx  > 

2030  Y 1 =Vgr i d  < 1 > 

2040  Y2=Ygr i d<Ny) 

2050  GINIT  200 .  ''260.  !  VERTICAL  PAPER 

2060  PLOTTER  IS  SOS/'HPGL" 

2070  PRINTER  IS  505 

2080  PRINT  " VS2 " 

2090  LIMIT  PLOTTER  505 , 0 . , 200 . , 0 . , 260 .  !  1  GDU  =  2  mm 

2100  VIEWPORT  22. ,85. , 19. , 122. 

2110  WINDOW  X 1 , X2 , Y 1 , Y2 

2120  FOR  1=1  TO  Nx 

2130  MOVE  Xgr i d< I > , Y 1 

2140  DRAW  Xgr 1 d< 1 > , Y2 

2150  NEXT  I 

2160  FOR  1=1  TO  Ny 

2170  MOVE  XI , Ygri d< I > 

2188  DRAW  X2 , Ygr i d  < I ) 

2190  NEXT  I 

2200  LDIR  0 

2210  CSI2E  2.3, .5 

2220  LORG  5 

2230  Y=Y1-<Y2-Y1 >*. 02 

2240  FOR  1=1  TO  Lx 

2250  MOVE  Xcoor d< I ) , Y 

2260  LABEL  X1abel*<I) 

2270  NEXT  I 

2280  CSIZE  3. , .5 

2290  MOVE  .5*<X1+X2),Y1-.06*<Y2-Y1> 

2300  LABEL  A* 

2310  MOVE  .5*<X1+X2),Y1-. 1*<Y2-Y1> 

2320  LABEL  C$ 

2330  CSIZE  2.3, .5 

2340  LORG  8 

2350  X=X1-(X2-X1 )*. 01 

2360  FOR  1=1  TO  Ly 

2370  MOVE  X, Ycoord< I ) 

2380  LABEL  Ylabel*<I> 

2390  NEXT  I 

2400  LDIR  PI/2. 

2410  CSIZE  3. ,.5 

2420  LORG  5 

2430  MOVE  X1-. 15*<X2-X1 >, . 5*<Y1+Y2> 

2440  LABEL  B* 

2450  PENUP 
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2460  PLOT  Pf <*> , Pdl <*> 

2470  PENUP 

2480  PLOT  Pf <*),Pd2<*> 

2490  PENUP 

2500  PLOT  Pf  <  *  >  , Pd3<  *  )  • 

2510  PENUP 

2520  PLOT  Pf <*>,Pd4<*> 

2530  PENUP 

2540  PLOT  Pf  <  * ) , Pd5<  * ) 

2550  PENUP 

2560  PLOT  Pf <*> , Pd6<*> 

2570  PENUP 

2580  PLOT  Pf<*),Pd7<*> 

2590  PENUP 

2600  PLOT  Pf <*>,Pd8<*> 

2610  PENUP 

2620  PLOT  Pf<*>,Pd9<*> 

2630  PENUP 

2640  PLOT  Pf <*),Pdl0<*> 

2650  PENUP 

2660  PLOT  Pf <*> , Pdl 1 <*> 

2670  PENUP 

2680  PLOT  Pf <*) , Pdl2<*> 

2690  PENUP 

2700  PLOT  Pf <*>,Pdl3<*> 

2710  PENUP 

2720  PLOT  Pf <  *  )  , Pd  1 4  <  * ) 

2730  PENUP 

2740  PLOT  Pf <*) , Pdl5<*) 

2750  PENUP 

2760  FLOT  Pf <*),Pdl6<*> 

2770  PENUP 

2780  PLOT  Pf<*>,Pdl7(*) 

2790  PENUP 

2800  PLOT  Pf <*) , Pdl8<*) 

2810  PENUP 

2820  PLOT  Pf <*>,Pdl9<*> 

2830  PENUP 

2840  PLOT  Pf <*>,Pd20<*> 

2850  PENUP 

2860  BEEP  500,2 

2870  PAUSE 

2880  PRINTER  IS  CRT 

2890  PLOTTER  505  IS  TERMINATED 

2900  SUBEND 
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